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We present several results obtained using the BlackHat next-to-leading order QCD program 
library, in conjunction with SHERPA. In particular, we present distributions for vector boson plus 
1,2,3-jet production at the Tevatron and at the asymptotic running energy of the Large Hadron 
Collider, including new Z+ 3-jet distributions. The Z+ 2-jet predictions for the second-jet Pj 
distribution are compared to CDF data. We present the jet-emission probability at NLO in W + 2- 
jet events at the LHC, where the tagging jets are taken to be the ones furthest apart in pseudora- 
pidity. We analyze further the large left-handed W polarization, identified in our previous study, 
for W bosons produced at high Pj at the LHC. 
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1. Introduction 

The dawn of the Large Hadron Collider (LHC) era brings renewed incentive to continue im- 
proving theoretical predictions of Standard-Model backgrounds to new physics searches. For many 
searches, including some channels for the Higgs boson and for dark matter particles, the signals 
will be excesses in jet + lepton or jet + missing Ej distributions. Such signals can be mimicked by 
Standard-Model processes; accordingly, a thorough and quantitatively reliable theoretical predic- 
tion is needed. This requires a calculation through next-to-leading order (NLO) in QCD. 

Leading-order (LO) computations, while an important first step, suffer from a strong depen- 
dence on the unphysical renormalization and factorization scales. At this order, they enter only 
through the strong coupling a s and parton distribution functions, uncompensated by any behav- 
ior of the short-distance partonic matrix elements. Because the QCD coupling is large and runs 
quickly, the absolute normalization of cross sections has a substantial dependence on scales. For 
reasonable scale variations, the dependence is of the order of ±40% for the V + 3-jet processes we 
shall study, with V a heavy electroweak vector boson. The dependence also grows substantially 
with increasing number of jets. At NLO, the virtual corrections introduce a compensating depen- 
dence on the scales. The scale dependence shrinks to ±10%, and we obtain a quantitatively reliable 
answer. Shapes of distributions can also show a dramatic scale dependence with poor scale choices. 
Some shapes do display noticeable "genuine" NLO corrections, independent of scale issues. 

NLO predictions for V + «-jet production at hadron colliders require several ingredients: 

• tree-level V + (n + 2)-parton matrix elements, which provide the LO contribution; 

• interference of one-loop and tree amplitudes for V + (n + 2) partons (virtual contribution); 

• tree-level V + (n + 3) -parton matrix elements (real-emission contribution); 

• a subtraction approximation capturing the singular behavior of the real-emission term; 

• the integral of the approximation over the singular phase space (real-subtraction term). 

These contributions must be convoluted with parton distribution functions, obtained from NLO fits, 
and integrated over the final phase space, incorporating appropriate experimental cuts. 
Schematically, we combine the contributions as follows, 



dObs 



■ J dx h2 M 2 J d<S> n Sobs C^ v+n + J d<S> n Sobs {O^ln + 
+ Jd® n+ i Sobs {G2™V+n+l - V^V+n+l) 



(1.1) 



where d<$> n denotes the V ±«-parton phase space; dx\p,f\fi the integral over the appropriate parton 
distributions, a sum over types being implicit; 5obs, the binning function for the desired distribution; 
c tree , the tree-level squared matrix elements; a 1 " 100 ?, the virtual corrections; a app , the approxima- 
tion to the real-emission contribution; and c^ app , the approximation's integral over singular phase 
space. The set of subtraction terms ensures that each of the terms in this equation is separately 
finite, and thus may be computed numerically. 
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We use the BlackHat program library [jl|, ||, ||, to compute the virtual corrections a 1 "' 0011 , 
and the SHERPA package [§] to compute a tree and the required approximation (a app and a^ app ). 
The approximation uses the Catani-Seymour dipole approach The phase-space integration is 
performed with SHERPA, implementing a multi-channel approach 

The BlackHat library implements on-shell methods for one-loop amplitudes numerically. 
Such amplitudes can be written as a sum of cut terms C n , containing branch cuts in kinematic 
invariants, and rational terms R n , free of branch cuts, 

A n = C n + R n . (1.2) 

All the branch cuts appear in the form of logarithms and dilogarithms, and can be written as a sum 
over a basis of scalar integrals — bubbles I' 2 , triangles I' 3 , and boxes l\, 

C n = I>4 + L C ' r 3 + I> 7 2- (1-3) 
(' (' i 

(Massive particles in the loop also require tadpole integrals.) We take all external momenta to be 
four dimensional, expressible in terms of spinors. The coefficients of these integrals, &i,c,-, and 
dt, as well as the rational remainder R n , are rational functions of spinor variables (in the form of 
spinor products). The BlackHat library computes these coefficients numerically, leveraging off 
recent analytic progress. In particular, it exploits generalized unitarity [pi 0]. We use Forde's ap- 
proach JTQ[ ] to compute bj and c,, making use also of the subtraction approach to integral reduction 



first introduced by Ossola, Papadopoulos and Pittau [11]. To obtain the rational terms we have 



implemented both loop-level on-shell recursion [12], and a "massive continuation" approach due 



to Badger Q13Q, which is related to the D-dimensional generalized unitarity Q14Q approach of Giele, 



Kunszt and Melnikov [15] 



One-loop matrix element computations can suffer from numerical instabilities. In Black- 
Hat, this problem is solved by detecting pieces of the amplitude which do not have a sufficient 
accuracy and recomputing them with higher precision using the multiprecision package QD [jig]. 
This approach has the advantage of solving the problem using the same approach for well-behaved 
points and for numerically unstable ones. As discussed in refs. [jl], |J, with a series of tests — the 
simplest of which checks whether the infrared divergences have the proper values — there is no 
need for a priori knowledge of what set of circumstances can lead to instabilities. In each con- 
tribution where precision loss is detected, BlackHat automatically switches to higher precision, 
regardless of the underlying cause. With on-shell methods this happens infrequently and therefore 
has only a mild effect on the overall computation time. 

We have previously used these software tools to provide the first phenomenologically useful 
NLO study of the production of a W boson in association with up to three jets [|3|,^j]. In this Contri- 
bution, we extend our previous studies with a more detailed look at the question of scale choices; 
at aspects of the polarization of Ws produced at high Pj\ and at a new distribution displaying the 
probability of emitting a jet into a rapidity gap. We also present the first NLO results on Z + 3-jet 
production at hadron colliders, in a leading-color approximation designed to be accurate within 
a few percent. In all cases, we decay the vector boson to leptons, W + — > Z + V/, W~ — > l~Vi, and 
Z — > l + l~ , using the appropriate vector boson linewidth. We include the virtual photon contribu- 
tion to l + l~ production. Other recent state-of-the-art NLO results may be found in ref. [17]. The 
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production of W + 3 jets has also been computed at NLO using a leading-color approximation and 
extrapolation [18, 19]. 



2. Scale Choices 

The renormalization and factorization scales are not physical scales. Physical quantities should 
be independent of them. A dependence on them is nonetheless present in theoretical predictions 
that are truncated at a fixed order in perturbation theory. At leading order, the dependence arises 
solely through a s and the parton distributions, respectively. We adopt the usual practice and choose 
the two to be equal, /i# = pLp = pi. NLO results greatly reduce the dependence compared to LO, but 
of course they do not eliminate it completely. We still need to choose this scale. We should expect 
a good choice for p, to be near a typical energy scale for the observable we are computing, in order 
to minimize the uncomputed logarithms in higher-order terms. However, multi-jet processes such 
as V + 2, 3-jet production have many intrinsic scales, and it is not clear a priori how to distill them 
into a single number. For any given point in the fully-differential cross section, there is a range of 
scales one could plausibly choose. For example, one might choose the same fixed scale p for all 
events. However, because there can be a large dynamic range in momentum scales (particularly at 
the LHC, where jet transverse energies well above Mw are common), it is natural to pick the scale 
p dynamically, event by event, as a function of the event's kinematics. 
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Figure 1: LO and NLO predictions for the second jet Et distribution in W - 
The only difference between the left and right panels is the scale choice: fi 



3 jet production at the LHC. 
= on the left and [i = Ht 
on the right. The former choice is clearly problematic and should not be used in phenomenological studies. 
The bottom panels show the LO and NLO predictions, varied by a factor of two around the central scale, and 
divided by the NLO value at the central scale. 



Previous studies (see e.g. refs. |2C, 21]) have used the transverse energy of the vector boson, 



Ej, as the scale choice. For many distributions at the Tevatron, this is satisfactory. With the larger 
dynamic range at the LHC, the choice becomes problematic. Indeed, for some observables, such as 
the transverse-energy distribution of the second-hardest jet in W + 3-jet production, shown in the 
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Figure 2: The NLO Pj distribution of the third jet in Z + 3-jet production at the Tevatron. For the left 
panel the scale choice jj, = Ej is used, and for the right panel ji = Hj/2. Although the two NLO results are 
compatible, the LO results have large shape differences, illustrating that jj. = Hj /2 is a better choice than 
jj. = Ej at the Tevatron as well. The lepton and jet cuts match the CDF ones ||2C|]. 

left panel of fig. |IJ, it goes disastrously wrong, leading to negative values of the distribution for Ej 
beyond 475 GeV. Even at the Tevatron, the scale choice \i = Ej is not necessarily a good one; for 
example, with this choice, the left panel of fig. || displays a large change in shape between LO and 
NLO in the Pj distribution of the third hardest jet in Z + 3-jet production. This difficulty reflects 
the emergence of a large logarithm ln(/i /E), where E is a typical energy scale, spoiling the validity 
of the perturbative expansion. 





(a) 

Figure 3: Two distinct W + 3 jet configurations with rather different values for the W transverse energy. 
In configuration (a) an energetic W balances the energy of the jets, while in (b) the W is relatively soft. 
Configuration (b) generally dominates over (a) when the jet transverse energies get large. 

To understand the problem with the scale choice /I = Ej, consider the two configurations 
depicted in fig. ||. In configuration (a), the W has a transverse energy larger than that of the jets, 
and accordingly sets the scale for the process. In configuration (b), the two leading jets roughly 
balance in Ej, while the W has much lower transverse energy. Here, the W scale is too low, and not 
characteristic of the process. In the tail of the distribution, we expect configuration (b) to dominate, 
because it results in a larger second-jet Ej for fixed center-of-mass partonic energy; contributions 
from higher center-of-mass energies will be suppressed by the fall-off of the partem distributions. 

Can we choose a scale that treats the different final-state objects more democratically? The 
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total partonic transverse energy, 

H T = £ 4 +E e T + £ T , (2.1) 

partons (' 

or a fixed fraction of it, is such a choice. As we can see in the right panels of figs. [T] and 0, this 
choice results in stable and sensible NLO predictions — and also in a relatively flat ratio of the 
NLO and LO predictions. For LO predictions, it is better to use such a scale when NLO results 
are unavailable. A similar type of scale choice, based on the combined invariant mass of the jets, 
has been motivated by soft-collinear effective theory [|22|]. Local scales associated with "branching 
histories" as used in parton showers have recently been studied for W + 3-jet production at LO [JH]]. 



3. Z + Jets at the Tevatron 

At hadron colliders, Z boson production manifests itself primarily in either charged-lepton 
pair production, or the production of missing transverse energy (when the Z decays to neutrinos). 
The latter process is an important background to a wide variety of supersymmetry searches (when 
no charged lepton is required), and to dark matter searches more generally. The l + l~ mode has 
a significantly lower rate, but it is an excellent calibration process, as the Z can be reconstructed 
precisely. It is also an excellent process for confronting NLO predictions with experimental data. 

We have computed the NLO Z+ 1,2, 3-jet production cross sections for the Tevatron (pp 
collisions at y/s = 1.96 TeV), with the Z decaying into a charged lepton pair. We applied the same 
cuts used by the CDF collaboration [p0|] in their measurement of these processes for Z — > e + e~, 



4 et > 30 GeV , E e T > 25 GeV , A/? e _ jet > 0.7 , 66 < M e+e - < 1 16 GeV , 
|T] jet |<2.1, \r) ei \<\, W 2 \<1 or 1.2 < \r] e2 \ < 2.8, (3.1) 

where the electron cuts apply to both electrons and positrons, and the jet cuts apply to all jets. We 
cut on the jet pseudo-rapidity r\ rather than CDF's cut on rapidity y; the two cuts coincide at LO 



but differ slightly at NLO. We employed three different infrared-safe jet algorithms [|23p , SISCone 
(with merging parameter / = 0.75), kj and natx-kj, all with R = 0.7. Production of an l + l~ pair 
can also be mediated by a virtual photon; we include these contributions as well, although they are 
suppressed by the cut on the lepton-pair invariant mass M e + e -. 

Fig. U shows how the Z+ 1,2, 3-jet cross section depends on a fixed scale }i, independent of 



the event kinematics, for the anti-^r algorithm and with the cuts (3.1). Here choosing /J. rj Mz is 
appropriate, because the cross section is dominated by low-/V jets. The upper three panels show the 
scale dependence of the cross section at NLO, compared to that at LO, in Z + 1-, Z + 2-, and Z + 3- 
jet production, respectively. They illustrate the lessened dependence at NLO. The bottom panel 
shows the ratio of NLO to LO results for all three cases, demonstrating the increasing sensitivity 
to scale variations at LO with increasing number of jets. This is expected, because there is an 
additional power of a s (n) multiplying the LO cross section for each additional jet. Accordingly, 
the impact of an NLO calculation also grows with the number of jets. The results for the kj and 
SISCone algorithms (not shown) are similar. 

Fig. H compares the theoretical predictions for the second-jet Pj distribution in Z + 2-jet pro- 
duction with data from CDF p0|]. CDF used the midpoint algorithm [^]. This algorithm is 
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Figure 4: The scale dependence of the cross section for Z+ 1,2,3-jet production at the Tevatron, for the 
anti-^7 jet algorithm using a leading-color approximation with rtf terms, as a function of the common renor- 
malization and factorization scale ji, with /Xo = Mz- The bottom panel shows the K factors, or ratios between 
NLO and LO results, for the three cases. 
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Figure 5: The second-jet Pj distribution for Z + 2 jets at LO and NLO compared against CDF data [E 



infrared unsafe for Z + 3-jets at NLO, so we use infrared-safe ones instead. Fig. |5| shows results 
for the anti-^7 algorithm; the other two algorithms yield similar results. It is worth noting that 
CDF did not attempt to "deconvolve" the hadronization corrections (estimated using Pythia) from 
their measured data; rather, they provided a table of hadronization corrections. This is helpful be- 
cause it will allow for future improvements to hadronization models to be taken into account in 
theoretical predictions. Accordingly, we have used these hadronization corrections to generate a 
complete prediction from the LO and NLO perturbative predictions. The hadronization corrections 
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are significant for low Pj, on the order of 20 % at 30 GeV, and become rather small at larger jet 
transverse momenta. As expected, the LO scale-dependence band is much larger than the NLO 
one. Excepting perhaps the last bin, the agreement between the NLO prediction and the data is 
quite good, especially given the different jet algorithms. 

Fig. ^ gives our predictions for the three jet Pj distributions in Z + 3-jet production, using the 
anti-^r jet algorithm. With the choice of scale /I = Hj/2, only minor shape changes are visible 
between LO and NLO, for all three distributions. The NLO plots are based on a leading-color 
approximation along the lines of refs. [^, ffl], except that pieces proportional to the number of light 
quark flavors (n f) are included. We expect this approximation to be valid to a few percent. 




First Jet P [ GcV ] 



Second Jet P. r [GeV] 



200 50 100 150 200 
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Figure 6: The LO and NLO Pj distributions for Z+ 3-jet production for the leading, second and third jet, 
for the anti-ftr algorithm and scale choice fi = Ht/2. The thin vertical bars in the top panels indicate the 
integration errors. 



4. W Polarization at the LHC 




Charged Lepton [ GeV | Neutrino E T | GeV ] 

Figure 7: The left panel shows the ratio of the charged-lepton Et distributions at the LHC for W + and W~ 
production in association with at least three jets, computed at NLO. The right panel shows the corresponding 
ratio for the neutrino Ej, or equivalently Ej. 

As noted in ref. [Q], at the LHC the Ej distributions of the daughter leptons show a surprisingly 
strong shape dependence on whether they come from a W + oraff", independent of the number 
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of jets. Fig. [7] shows the ratio of the NLO transverse energy distributions for the W ± boson decay 
products in inclusive W + 3-jet production at the LHC, charged leptons in the left panel and neutri- 
nos in the right panel. The differences between W + and W ~~ distributions are quite dramatic. The 
left panel shows a large ratio for W + to W~ at small Ef which declines at larger Ej. In contrast, the 
corresponding ratio for the Ej, or equivalently the missing transverse energy Ej in the event, starts 
somewhat smaller but increases rapidly with Ej. The significant difference in behavior between 
W + and W~ suggests a means for separating W bosons produced in top quark decays from those 
produced from light quarks; the Ws from top decays do not exhibit a similar phenomenon. 
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Figure 8: The LO and NLO predictions for polarization fractions of the left-handed, fi (top curve), right- 
handed fn (middle curve) and longitudinal /o (bottom curve) fractions for W + 2 jets at the LHC. The left 
panel gives the polarization for W + and the right panel for W~ . For high transverse momentum, Pjyv, the 
W bosons become predominantly left-handed. 

This disparate behavior is explained by a net left-handed polarization for both W + and W~ at 
high transverse momentum. This effect is easily visible at LO, and it does not get washed out at 
NLO. In fig. II we give the fraction of W bosons in each of the three polarization states, left-handed, 
right-handed and longitudinal (/l,/r,/o> respectively) for W + 2-jet production at the LHC, at both 
LO and NLO. As seen in the figure, at high transverse momentum the W ± bosons are preferentially 
left handed. Although the cross-sections for W + and W~ are rather different, their polarizations 
are nearly identical. Interestingly, we also find that when the Ws have a transverse momentum of 
more than 50 GeV, the polarization is quite independent of the jet transverse energy cuts. With 
W^ 1 bosons left-hand polarized at large E^, the W + tends to emit the left-handed neutrino forward 
relative to its direction of motion (resulting in a larger transverse energy) and the right-handed 
positron backward (smaller transverse energy). In contrast, the W~ prefers to emit the left-handed 
electron forward. At high Ej, such decays produce an enhancement in the neutrino Ej distribution 
and a depletion in the charged-lepton distribution, for W + relative to W~, consistent with the results 
displayed in fig. [7[ We note that this phenomenon is distinct from the well-known dilution of the 
W rapidity asymmetry at the Tevatron, when passing to the decay lepton, which can be explained 



using angular momentum conservation solely along the beam axis [|25|]. 

5. Emission into Rapidity Gaps 

In previous work [pT|], we provided the first NLO study of the probability of emitting a third jet 
in W + 2-jet events, as a function of the rapidity interval between two leading-£y jets at the LHC. 
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This distribution was studied earlier at LO at the Tevatron and compared to CDF data [26]. Jet 



emission probabilities are relevant to Higgs searches in vector-boson fusion Q27|], in which color- 
singlet exchange leads to a paucity of jet radiation in the central region between two forward tag 
jets. On the other hand, QCD backgrounds with color exchange, as in W + 2-jet production, will 
generally lead to significant jet radiation. 



LO: (W +3jets + X) / (W +2jets + X ) 
NLO: ( W" + 3 jets + X ) / ( W" + 2 jets + X ) 




Figure 9: The ratio of the inclusive W~ + 3-jet cross section to the W~ + 2-jet cross section as a function of 
the pseudorapidity separation At] between the two most widely separated jets that pass the cuts. The solid 
(black) line gives the NLO result, while the dashed (blue) line gives the LO results. 



To mimic vector-boson fusion searches, however, the appropriate tag jets are not the two hard- 
est ones (by Ej), but rather the two most separated in pseudorapidity. Therefore, in fig. ^ we 
present the ratio of the W~ + 3-jet cross section to the W~ + 2-jet cross section as a function of 
the pseudorapidity separation At] between the two most separated jets. The emission probability 
rises roughly linearly with At] . The NLO result is somewhat less than the LO one at large At] . 
(The ratio for W + is quite similar.) This plot is similar to one for Higgs production in association 



with jets [28], obtained from high-energy factorization considerations. It would be interesting to 
compare results obtained in this way to NLO results for the same quantities. 



6. Conclusions 

In this Contribution we presented some new results for W + 3-jet production obtained from 
BlackHat combined with SHERPA, expanding on earlier scale-dependence studies [||, We 
also demonstrated that W bosons produced at large Pj are indeed polarized left-handed, explaining 
an asymmetry between W + and W~ in the transverse energy distributions of the daughter leptons. 
Because Ws from top decays do not exhibit this polarization effect, it may prove effective for 
distinguishing such Ws from ones produced by light quarks. We presented the first NLO study 
of the probability of emitting a third jet between the two most widely separated jets in W + 2- 
jet production. We also presented the first NLO results for Z+ 3-jet production. We observed that 
even at the Tevatron, choosing the renormalization and factorization scale to equal the vector boson 
transverse energy is not a particularly good choice, as it induces large shape changes between LO 
and NLO. 

A publicly available version of BlackHat is in preparation and is currently being tested in 
diverse projects (see e.g. ref. p9|]). This version uses the proposed Les Houches interface for one- 
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loop matrix elements. It has been tested with both C++ and Fortran clients. The public version will 
provide all processes that have been carefully tested with the full BlackHat code. 

In the more distant future, the next benchmark process for BlackHat + SHERPA is the pro- 
duction of a W boson in association with four jets at NLO. Using the techniques described above, 
the virtual part of the NLO cross section seems within reach. Computing the real emission ma- 
trix elements, and integrating them over the seven-particle phase space (including the decay of the 
vector boson) appears to be rather challenging with the current tools, due to the large number of 
integration channels. It is interesting to note that in this case the bottleneck no longer seems to be 
the virtual contributions to the cross section. 

The results summarized here are indicative of the type of physics that can be carried out using 
BlackHat in conjunction with SHERPA. We look forward to comparing predictions from these 
tools to the forthcoming LHC data. 
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